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Abstract 

A basic concept to calculate physical features of non-ideal plasmas, such as optical prop- 
erties, is the spectral function which is linked to the self-energy. We calculate the spectral 
function for a non-relativistic hydrogen plasma in GW-approximation. In order to go beyond 
GW approximation, we include self-energy and vertex correction to the polarization function 
in lowest order. Partial compensation is observed. The relation of our approach to GW and 
GWT calculations in other fields, such as the band-structure calculations in semiconductor 
physics, is discussed. From the spectral function we derive the absorption coefficient due to 
inverse bremsstrahlung via the polarization function. As a result, a significant reduction of 
the absorption as compared to the Bethe-Heitler formula for bremsstrahlung is obtained. 

1 Introduction 

Spectroscopy can serve as a versatile tool to infer properties of dense and hot plasmas from the 
emitted radiation^ |2]. In non- ideal plasmas, where the coupling parameter can exceed unity, 
properties often deviate significantly from their form in ideal plasmas due to the importance of 
interaction effects With the advent of femtosecond laser pulses 0], it has become possible 
to produce non-ideal plasmas with table-top systems. In particular, it is nowadays possible to 
create conditions similar to those in the center of astrophysical objects such as the sun or giant 
planets. Effects such as dynamical screening, dissolution of bound states, Pauli-blocking, and the 
importance of collisions have been observed : 3: . Taking proper account of interaction effects is a 
challenge to any theoretical description of non-ideal plasmas. 

Many-body perturbation theory presents a toolbox to determine various properties of non-ideal 
plasmas [51|H1[7]- The use of Green's function techniques allows for a systematic and intuitive 
consideration of many-body effects. In particular, one-particle properties can be obtained from 
the one-particle Green's function or its spectral representation, the one-particle spectral function. 
A number of important mechanisms in dense plasmas such as dynamical screening of the Coulomb 
interaction can be described by partial summation of certain diagrams. A particularly successful 
concept is the quasi-particle picture 00 [7|- However, with increasing coupling the quasi-particle 
pictures breaks down. This shows up in a broadening of the spectral function. To go beyond 
the quasi-particle picture in a consistent way poses serious problems of self-consistency. Vertex 
and self-energy corrections have to be taken into account on the same footing to obey sum rules 
and other exact known properties. Using Ward identities |S] or the Kadanoff-Baym scheme of 
the conserving vertex [5] enables one to construct consistent sets of diagrams, alas the resulting 
integral-equations are complicated to solve. 

The GW approximation is a particularly scheme for the self-energy approximating it by a 
product of the Green's function G and an effective interaction W ^H]- Being introduced by Hcdin 
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in 1965 m, it has a long history ol applications, which is reviewed in Ref. ^2] and ^S]- However, 
the full self-consistency implied in Hedin's original proposal was not carried out so far. Instead, the 
ideal Green's function or a quasi-particle picture was used. In this manner, the GW approximation 
was successfully applied to determine shifts in the ionization potential, approximations for the 
equation of state, and an effective interaction to describe bound states in a medium |7j. Also, the 
band structure of different types of materials, like semiconductors [E], alkali an d transition 
metals jlfi) was determined. Due to the efficiency of modern computer technology, it has recently 
become feasible to address the self-consistency implied in Hedin's scheme to some extent, see 
e.g. |171 1181 ITUl I20| . Some of these calculations improve solely on self-energy corrections, while 
others take into account vertex contributions as well. It is customary to call the later GWT 
approximations. In this paper, we perform calculations for the optical properties of a hydrogen 
plasma at solar core conditions, as an important example for astrophysical plasmas. Typical 
parameters of the solar core plasma are temperatures of about T = 100 Ry ~ 1360eV/&B and 
particle densities reaching n ~ 7T0 24 cm~ 3 . These parameters justify the model of a classical (non- 
degenerate) and weakly coupled plasma, being characterized by values of 9 — k^T / ~ 10 for 
the degeneracy parameter and V = (Ze 2 /47reo)/(3/47rn) 1 / 3 fcBT ~ 0.03 for the coupling parameter. 

Central to the description of optical properties of plasmas is the dielectric function e(A;,w) 
PQ. Within the frame of the approach presented here, it can be obtained from the polarization 
function, which is a member of Hedin's set of equations. In an earlier paper |21|. the suppression 
of the bremsstrahlung cross section due to successive scattering has been treated by a taking 
into account self-energy and vertex corrections. However, these corrections were only considered 
within a one-loop approximation. In this way, dynamical screening effects were neglected. It is the 
objective of this communication to present results with an improved one-particle spectral function 
as obtained with the GWT approximation. 

In Sec. [3 we will review the GW approximation making use of Hedin's equations. Sec. 0] 
presents an illustrative example for the GW° approximation. Implications for the absorption 
coefficient are discussed in Sec. El A discussion and conclusions will be given in Sec. 

If not otherwise indicated, we apply the Rydberg system of units where h = fee = 1 , e 2 = 
2 , e = 1 /4-7T , and m e = 1/2. 



A convenient starting point of our approach are Hedin's equations 22 j. It is a closed set of 
equations relating the full Green's function G, the non- interacting Green's function Go, the self- 
energy S, the dynamically screened interaction W a b I10| - the polarization function n, and the 
vertex function T. In detail, the full Green's function is given by 



Here and in the following, the shorthand notation (1) = (ri, ti, <Ti . . .) for spatial variables r, 
imaginary times r as well as quantum numbers such as spin is used. The self-energy S Q (12) is 
obtained from the dynamically screened interaction and the vertex function according to 



2 The GWT approximation 




(i) 





The dynamical screened interaction is given via 





by the polarization function 




(4) 
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Finally, the vertex function obeys a Bethe-Salpeter like equation 

r„(12,3) = 5(12) 5(13) + J2 y^(4567) ^|^ G b (46)G b (75)r b (67,3) . (5) 

The set of equations represents a perturbation expansion in terms of the screened interaction 
W a b and is expected to show better convergence properties compared to an expansion in the 
inter-particle interaction V a b- 

By truncating the set of equations Q-® on a certain level, various approximations for these 
quantities can be defined. A particular simple approximation is the GqW° approximation |11II22| . 
It is obtained by taking the bare vertex T(12, 3) = 5(12)5(13) and inserting non-interacting Green's 
functions into the expression for the polarization function. This leads to a self-energy given by 
Y^ w (12) — iG a (12)W aa (21). Next, the Green's functions in this expression as well as in the 
screened interaction W are taken as the free Green's function Go- This expression turns out to 
be quite successful. However, in dense plasmas as well as in a number of materials in solid state 
physics, there is need to improve beyond this simple approximation. 

Such an improvement is the GW° approximation which was studied in Ref. [231 f° r the electron 
gas at T = and in Ref. for the solar core plasma. In the latter case, a considerably broadened 
quasiparticle and a featureless behaviour at the plasma frequency was found. No plasmon-like 
satellite structures survived the partial self-consistency treatment. 

Another straightforward extension is a self-consistent solution of Dyson's equation and the 
screened interaction while keeping the vertex function in lowest order |18| . This will be termed 
GW approximation in this paper. In Ref. |18| . an increase in the quasiparticle bandwidth and a 
featureless satellite structure is found in contradiction to experimental evidence. Also, such an 
approximation leads to a drastic violation of the f-sum for the inverse dielectric function. The 
results indicate the importance of vertex corrections [2SI - Calculations for a one-dimensional 
semiconductor lead to similar conclusions. Large cancellations between self-energy and vertex 
corrections have also been found in Ref. |27l 1281 1291 1301 l3*T] . 

Recently, a number of approximation schemes to take into account vertex corrections have been 
proposed $X ■ I n this work we will apply a sequence of approximations described in the following: 

Taking the self-energy in the GVK-approximation, the functional derivative occurring in the 
vertex equation JSJ) yields in lowest order of W aa 

ggj = ^(12)5(14)5(25)5^ . (6) 

Note that the screened interaction is taken in the one-loop approximation for II, the so-called 
random-phase approximation (RPA) . This term leads to a ladder approximation for the vertex in 
terms of the dynamically screened interaction 

T a (12,3) = 5(12)5(13) + /rf(67)IF a ° a (12)G Q (16)G a (72)r a (67,3) . (7) 



Already this equation is challenging to solve |33|. even in Shindo approximation 34 . As a result, 
the improved self-energy in this approximation is given in second order of the dynamically screened 
interaction by the term studied in Ref. |17| . 



(12) = G (12)W°(21) + /d(34)G (13)Wj (23)G (34)TC(41)G (42) . (8) 



Also, the improved polarization function is given besides the loop diagram by an exchange diagram 
with respect to W° 



l£Ta (12) = G„(12)G a (21) + /d(34)G (13)G a (41)W o u o (43)G a (24)G a (32) . (9) 

At this stage the first iteration of Hedin's equations GJ-10 I s completed. We will not go beyond 
this first iteration, in particular the screened interaction potential W aa is kept on the level of RPA. 
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To address the deviations from the quasi-particle picture, the introduction of the one-particle 
spectral function A a (p,u>) is convenient. With its help, the spectral representation of the full 
Green's function is given by 

r, i a [°° duj A a (p,uj) 

Ga(p,z) = / , (10) 

J.o^Tr z-w 

p and z denote momentum and energy (frequency) of the particle. The Green's function G a (12) 
is obtained from the Green's function in momentum-frequency representation G a (p, z) by means 
of Laplace transform. The quasi-particle approximation itself can be stated as a <5-likc form of the 
spectral function 

A a {p,Lu) = 2tt<5 (w - E a (p)) , (11) 
where the energy E a (p) is obtained as 

2 

E a (p) = +ReX a (p,E a (p)) . (12) 

Note, that the spectral function itself obeys a normalization relation 

/ — A a (p,oj) = 1 . (13) 
Furthermore, sum rules for the first and second moment of the spectral function are known |23) . 

&uuA a (jp,uj) = Ef F , (14) 



/OO pOO 
duu 2 A a {p,u>) = / dwImES(p,w) + (Eaip)™) 2 , (15) 
-oo J — oo 

where the index HF refers to the energy in Hartree-Fock approximation and the index c indicates 
the use of the correlated self-energy. From the spectral function, a number of thermodynamic 
quantities can be obtained, e.g. the one-particle density [7113*5]. 

/j 3 r 
(2tt)3 7 7^ fa(u)A a (p,Uj) , (16) 

with the distribution function / (w) of particles of species a. 

Optical properties, which are under consideration in this work, can be obtained from the 
well-known relation between the dielectric function e(q,ui) and the polarization function Eq. flljl. 
i.e. 

e(q,w) = 1 - J2 V aa(q)^aa(q,i0) , (17) 

a 

with the interaction potential V aa {q). In particular, the absorption coefficient, which gives the 
attenuation of electromagnetic radiation traversing the plasma is given by the imaginary part of 
the long wavelength limit of the dielectric function as 

a(u)) = -Imefg -> 0,tA (18) 
c 

relation that holds for wavelengths long against interatomic distances and frequencies high com- 
pared to the plasma frequency oj p i = 4(7rn) 1 / 2 . For details, we refer to Refs. [Tl I24| . 
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3 Applications of GW approximation for the self-energy 



The GW approximation has been very successfully used for a broad variety of problems in many- 
particle physics for a long time. Here we would like to mention only a few examples. 

One of the first applications in solid state physics was the calculation of band structures com- 
bining density functional theory (DFT) and the GW-method. Northrup et al. ^5] showed how 
this approach improves the band-gap problem of the local density approximation (LDA), which 
predicts too small band-gaps for most materials, whereas the use of GW helps to decrease the de- 
viation of the theoretical from the experimental value below 0.1 eV in the case of silicon, compared 
to about 1 eV as obtained from pure LDA-DFT. Similar results were obtained for many different 
semi-conductors and insulators |13U37) . Optical properties of highly excited semiconductors using 
a combined GW and T-Matrix approach have been studied extensively by Schmielau et al., see 
Ref's |38l Kffll 1471] . Also in the case of finite systems, the GW approximation has been applied suc- 
cessfully. For the case of the Na4 tetramer, good agreement between the theoretically predicted 
photoabsorption and experimental data was obtained |41j . 

Nuclear and quark matter have been investigated using techniques very similar to the approach 
presented here. The problem of dynamical chiral symmetry breaking, which cannot be described in 
a perturbative calculation of quark propagators, was demonstrated to emerge in a nonperturbative 
approach known in quantum field theory as rainbow-ladder approximation |42l I48j . 

Fehr et al. 03] and Wierling et al. [23] investigated spectral properties of electrons in nonideal 
plasmas making use of the GW approximation. Whereas Fehr applied the perturbative GoW° 
approach (c.f. section |5J) to both equation of state and optical properties, Wierling already used 
a self-consistent GW° approximation which is also used here in the following section. 



4 GW° approximation for the solar core plasma 

In the GW° approximation, the set of Hedin's equation reduces to a self-consistent solution of 
Dyson's equation with the correlated part of the self-energy, 

A c (p, UJ + IVj) = 2 2 ' (19) 

[uj — p 2 - Re S c (p, uj + i-q)] + [Im S c (p, u + i?])] 

oo oo 

Y, e (p,uj + tr)) = - / / — -tt— Ime RPA +i0) l+n B (w) — — - j ■ 



n 



(20) 



As indicated by W°, the dielectric function e(q,oj) is taken in random phase approximation (RPA) 
and the Bose function is given as ns(u) = (exp (—oj/T) — 1) . Here, as in the following, electrons 
labeled by a — e are considered. Note that the Hartree-Fock self-energy, which appears as an 
additional term in Eq. l]2()[l is small for the considered parameters and is henceforth neglected. 

We give results for a hydrogen plasma using the conditions at the solar core center. We ignore 
other ions such as helium etc. in these exploratory calculations. In Fig. ^ the self-consistent 
spectral function calculated in GW° approximation, i.e. the solution of the set of Eqs. I|19|l and 
(J20J, is shown for a fixed momentum of p = O^lag 1 . The grey curve is the initial ansatz for 
the spectral function, i.e. the input for the r.h.s. of Eq. 1201) . It has been chosen of Lorentzian 
form with a width of 7 = 10 Ry. The first iteration is given by the dashed-dotted curve. The 
peak of the spectral function is shifted to smaller frequencies and the function is asymmetrically 
broadened. The second and third iteration give only minor modifications to the first iteration, the 
forth iteration (not shown) does not vary significantly from the third iteration. 

For comparison, we also show the GqW° approximation, i.e. the first iteration of Eq. H20|l 
starting from a spectral function of vanishing width (delta-distribution) . In this case, there is no 
quasiparticle peak at the quasiparticle resonance wqp = p 2 , but four satellites and minima appear. 
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Figure 1: Electron spectral function in GW° for solar core conditions (n c = 7 ■ 10 24 cm~ 3 ,T c = 
100 Ry ~ 1360 eV) at momentum p — 0.21<Zg . Dashed-dotted curve: first iteration of GW°, 
(quasiparticle approximation [QPA]). Grey solid curve: Ansatz for the spectral function in the 
iterative solution of GW°. Dotted curve: 1st iteration, dashed: 2nd iteration, solid: 3rd iteration. 
The 4th iteration is not distinguishable from the 3rd iteration. 



The latter are due to sharp peaks in the response function Im e R p A (&, u>) and have been observed 
earlier gSlEBI- 

As a result of this calculation, we observe, that the self-consistent calculation leads to a spectral 
function that is physical easily understandable, i.e. it contains a broadened and shifted quasipar- 
ticle resonance. However, the signatures of collective effects, such as the dynamical screening, 
which are present in the GqW° approximation (dash-dotted curve in Fig.^l vanish completely in 
the self-consistent result. In contrast, the GqW° result does not contain a broadened quasiparticle 
peak but is completely determined by the behaviour of the response function, which contains the 
collective excitations of the plasma (plasmons), as also shown by Fehr et al., see Ref. |47j . 

Although convergence is already achieved after 3-4 iterations of Eqs. 119(1 and 1(20(1 . the calcula- 
tion is too time-consuming to be used in the computation of physical observables such as equation 
of state or the polarization function, the latter involving convolution integrals over two spectral 
functions in both momentum and frequency domain. Therefore, we make use of some further 
approximations. First, we replace the full RPA-like screened interaction by a one-loop approxima- 
tion, which takes into account the scattering among particles in Born approximation. In this case 
we only consider electron-ion scatterings. Ions are treated in adiabatic approximation and the 
interaction is mediated by a statically screened potential of Debye-type, i.e. V e i(q) = %ir/(q 2 + k 2 ), 
where k = (87T7J/T) 1 / 2 i s the inverse Debye screening length. This approximation is justified by 
the fact that the corresponding absorption cross-section leads to the nonrelativistic limit of the 
Bcthe-Heitle formula for inverse bremsstrahlung in the limit of vanishing k, as shown in |21j . 

For small k, the main contributions to the integral in Eq. ((20(1 come from small q. Therefore, 
we can neglect the shift in the momentum variable in the self-energy on the r.h.s. As shown in 
Ref. [21| . one then obtains a particularly simple equation for the self-energy in this approximation, 



G 




CO [Ry] CO [Ry] 

Figure 2: Self-energy (real and imaginary part), dispersion relation and spectral function calculated 
self-consistently by solving Eq. pifl (left) and quasiparticle approximation (right), i.e. the r.h.s of 
Eq. J2U evaluated with = 0. Parameters:n c = 7 • 10 24 c m - 3 ,T o = 100 Ry ~ 1360 eV. The 

momentum is fixed at p = 0.21 ag 1 . 



which we will denote by the suffix GW°s, 



k 2 +p 2 - lu - irj + S GW s (p,uj + ir]) - 2inJ uj + ir] - E GW ° s (p,w + irj) 



S Gw0s (p, W+l?? )--87r^ 

K 

(21) 

which can be solved by standard root-finding algorithms. Fig. 2 shows the solution of Eq. Q21|l for 
the same set of parameters as used in the RPA calculation (Fig. To the left, the self-energy 
(imaginary and real part), as well as the dispersion relation, and the spectral function as obtained 
from the numerical solution of Eq. Q21[) are shown. The spectral function is centered around the 
quasiparticle energy. Its overall shape is given by the imaginary part of the self-energy, while the 
real part determines the dispersion. In the self-consistent solution, the dispersion relation contains 
only a single root near the quasiparticle energy, while there are three roots in the quasiparticle 
approximation (r.h.s of Fig. 2). The upper and lower root give rise to the two satellites in the 
corresponding spectral function, while the central root does not yield a resonance due to the large 
imaginary part of £ G ° W s around u = 1 Ry. The separation of the satellites from the quasiparticle 
energy is approximately given by the plasma frequency. This is a general feature, which has been 
observed in earlier calculations carried out at lower particle densities and temperatures |21| . 

In Fig. |21 we compare the spectral functions obtained from the two-component RPA and the 
simplified Born approximation containing only the scattering of electrons on fixed ions (Lorentz- 
Plasma). Since electron-electron collisions are neglected in the latter case, the corresponding 
spectral function is narrower than in the RPA calculation. Furthermore, in the one-loop ap- 
proximation there are no plasmon degrees of freedom which leed to further damping in the RPA 
calculation. On the other hand, the Lorentz plasma calculation already contains all characteristic 
features of the RPA spectral function. Namely the asymmetrical broadening and the shift of the 
peak towards lower energies show up. These features become less pronounced at higher momenta, 
where many-particle effects arc expected to be less important. With increasing momentum, the 
width of the spectral function decreases, while at the same time the height of the quasiparticle 
peak increases, as can be seen in Fig. 0] . This follows naturally from the normalization condition 
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Figure 3: Self-consistent spectral function computed in full RPA (two component plasma) (dashed 
curve) and in static Born approximation including only electron-ion collisions (solid curve). The 
neglectance of e — e collisions leads to a smaller width of the spectral function. Solar core param- 
eters: n = 7 ■ 10 24 cm" 3 , T = 100 Ry. 



for the spectral function Eq. (| 1 31) - In conclusion, we have shown that our simplified Lorentz plasma 
model for the self-energy leads physically intuitive one-particle spectral function which includes 
the electron-ion collisions in a consistent way. 

5 Applications of GWr-approximation for the self-energy 

In the last years, substantial effort has been made to go beyond the GVF-approximation and to 
include the vertex to some extent. In a pioneering paper, Takada |20j demonstrated the feasi- 
bility of self-consistent GWT approximation. He used an ansatz for the vertex-function which 
fulfils certain sum-rules and conservation laws. As a result, he showed the subtle cancellation of 
contributions from self-energy and vertex-corrections to spectral properties of charged particles, 
exemplified for the damping of plasmons in Al. Very recently, Ziesche |48| has reviewed the cal- 
culation of direct and exchange contributions (vertex correction) to the on-shell self-energy of the 
homogenous electron-gas. 

For plasmas, Vorberger et al. |35j systematically studied contributions to the equation of state 
beyond Montroll-Ward, including all exchange terms (vertex corrections). In the same spirit, we 
will now investigate the lowest order vertex correction to the self-energy in second order of the 
interaction, given by the following diagrammatic expression. 




The second diagram is the lowest order diagram which does not appear in the GW approximation, 



8 



whereas the first diagram is automatically included. The base lines are free electron propagators, 
while the upper loops are composed of ionic propagators taken in the adiabatic approximation. 
To ensure convergence, a statically screened potential V C i(q) is used. The numerical evaluation 
of Eq. 11' I'D is shown in Fig. [5J The dashed and dotted curves give the first and second iteration 
respectively of the GW°s approximation, the solid curve is the self-consistent result. For compar- 
ison, the first order vertex correction as obtained from the exchange diagram in Eq. 112211 is shown 
in the inset. The vertex term gives at most a 20% reduction of the self-energy. 

6 Modification of the absorption coefficient 

As an illustrative example, we calculate the impact of the broadened electron spectral function on 
the absorption coefficient due to inverse bremsstrahlung. The absorption coefficient for radiation 
in a plasma can be obtained from the dielectric function e(q,v), Eq. (fT%| . According to Eq. (@J), 
nonideality effects enter via self-energy corrections of the full Green's functions and by vertex 
corrections beyond the bare vertex. Within the GW approximation, only self-energy corrections 
are considered, while the GWT also accounts for vertex terms in the polarization function. 

6.1 Self-energy corrections 

Using the self-consistent spectral function obtained in simplified GW° (Eq. (|21|l ) to calculate the 
polarization function, leads us to the result plotted as solid curve in Fig.EI We normalize our result 
to Kramers' formula corrected by a Gaunt-factor in Born approximation |49U5l)| . which corresponds 
to the non-relativistic limit of the Bethe-Heitler cross section for inverse Bremsstrahlung In 
the infrared limit (to — » 0) the Born approximation shows a logarithmic divergence. 

The result for solar core parameters is shown in Fig.^J For small frequencies, a reduction of the 
free-free absorption of about 15% as compared to Born approximation is observed, while at high 
frequencies our approach converges to the Born approximation. At intermediate frequencies, we 
obtain a slight enhancement of our result relative to the Born approximation. These characteristics 
have already been observed in earlier calculations for lower densities and temperatures We 
compare our results to calculations which use parametrized spectral functions of Gaussian shape 
in the polarization loop (dashed and dotted curve). The width a is given. Also in this case, an 
enhancement is observed which decreases with decreasing width a as can be seen by comparing 
the dotted (cr = 6.5 Ry) and the dashed curve (a — 13 Ry). The value a = 6.5 Ry gives a Gaussian 
spectral function of similar shape as the self-consistent calculation at small momenta. As shown 
in Fig. 0J the self-consistently calculated spectral function converges to a quasi-particle resonance 
at large momenta, i.e. large frequencies. This behaviour leads to the faster convergence of the 
absorption coefficient to the Born result as compared to the calculation using Gaussians with 
frequency and momentum independent widths. Remember that the Born result for the absorption 
coefficient is obtained by inserting delta- like spectral functions in the polarization function Eq. . 

6.2 Vertex corrections 

We have shown that the account of a broadened one-particle spectral function leads to a suppres- 
sion of the infrared behaviour of the inverse bremsstrahlung spectrum. However, the calculations 
have been carried out on the level of GW°-approximation. Vertex corrections contribute to the 
polarization function in the same order with respect to the statically screened potential. More- 
over, these vertex corrections tend to cancel the self-energy correction to some extent. In order 
to account for both types of corrections on the same footing, we focus on contributions from the 
vertex equation, Eq. (JSJ, up to first order in the screened interaction in the polarization, Eq. Iprjl. 




(23) 
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In Fig. we compare the suppression of the absorption coefficient obtained from the inclu- 
sion of self-energy only (solid curve) and with included vertex correction in the polarization loop 
(dashed curve), c.f. Eq. J2SI- At the present parameters, the vertex correction is the most dom- 
inant term. Over the whole range of frequencies considered here, we obtain a suppression of the 
absorption coefficient as compared to the Born approximation. Also, the enhancement of absorp- 
tion at frequencies around to — 10 Ry, which was observed in the calculation using only self-energy 
corrections (solid curve), vanishes if the vertex is taken into account (dashed curve). This shows 
the importance of vertex corrections. 

Here, the vertex correction is only taken in lowest order of the interaction and density, while 
the self-energy based result contains a summation to all orders of density. In order to compare 
self-energy and vertex contributions to the modification of the absorption spectrum in a fully 
consistent way, one would have to go beyond the perturbative calculation presented here and solve 
for the vertex equation Eq. JSJ using at least free-particle propagators. This task goes beyond the 
scope of this paper, where only the lowest order corrections are to be studied. Finally, we remark, 
that in the mentioned earlier calculations presented in Ref. the vertex correction modifies 
the self-energy result only to a minor extent, i.e. the enhancement at intermediate frequencies is 
reduced by 40%. This is due to the lower density (10 19 cm' 3 ) used in that work. In all calculations, 
the high frequency behaviour converges nicely to the Born result. 

For the infrared part of the spectrum, the behaviour of the absorption coefficient needs further 
considerations, going beyond the present work. As shown in Fig. [TJ the self-energy correction in the 
polarization function reduces the absorption coefficient by a constant factor of about 15% in the 
low frequency limit, but does not regularize the infrared divergence of the Born approximation. 
The vertex term, on the other hand, induces a large suppression so that in the low frequency 
limit higher orders of the vertex correction have to be considered. Furthermore, it is well-known, 
that below the plasma frequency (7 Ry for the solar core parameters) dynamic screening plays an 
important role. This needs a further summation of diagrams going beyond the present scheme. The 
impact of dynamical screening on bremsstrahlung emission has been discussed by Ter-Mikaelyan 
in Ref. [S3- 

7 Conclusions 

A systematic treatment of optical properties in non-ideal plasmas is possible in the framework 
of Green's function methods. Corrections beyond the quasi-particle pictures can be generated 
using Hcdin's equations. However, a consistent solution of Hedin's equations is a formidable task. 
Here, we considered a simplified set of equations. Notably the GWo-appioximatioii, mainly used 
in solid state physics, has been shown to lead to sensible results also in the field of plasma physics. 
Results for the self-energy and the spectral function are presented. Plasmon-like structures, present 
in perturbative calculations as performed in Ref. |47| vanish completely. An asymmetrically 
broadened and shifted spectral function is obtained. Furthermore, it is shown, that a simplified 
model, where the dynamically screened interaction is approximated by static screening leads to 
similar results for both self-energy and spectral function as compared to the full RPA results. 
Deviations can be understood as a consequence of neglecting collisions among particles of equal 
species. 

The simplified GW° approximation for plasmas is now available for a broad range of parame- 
ters, i.e. density and temperature. Also, the correct asymptotic behaviour, i.e. convergence to a 
delta-like quasiparticlc resonance at large momenta, is obtained. In this paper we used parameters 
corresponding to the solar core, whereas Ref. |21| contains similar results for lower density and 
temperature. 

Vertex corrections to the self-energy have been studied on a perturbative level. So far, only 
the high frequency behaviour of the exchange diagram second order in the screened interaction 
was evaluated. It leads to 20% reduction of the self-energy with respect to the second order of 
GW°. 
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The availability of the one-particle spectral function for any momentum and frequency makes 
it possible to use it in calculations of further physical observables as the equation of state and 
the optical properties. Here we focused on the absorption of electromagnetic radiation due to 
free-free transitions, which is obtained from the polarization function. Insertion of the broadened 
particle propagators in the one-loop approximation for II leads to a significant modification of the 
absorption spectrum at low frequencies. Besides the suppression in the infrared, enhancement of 
absorption is observed at intermediate frequencies. 

Since the one-loop approximation using full propagators is an inconsistent approximation with 
regards to conservation laws, such as Ward identities, we investigated the vertex-correction inside 
the polarization loop in lowest order. It turned out, that the vertex correction is by far the 
most important correction to Born approximation, i.e. the self-energy effects contained in the 
full spectral functions, are dominated by the vertex correction. However, this is not a general 
feature as becomes clear if comparing to calculations for other parameters in Ref. 21 . On the 
other hand, a consistent comparison of both self-energy and vertex corrections on the same level of 
approximation necessitates the summation of vertex-corrections by solution of the Bethe-Salpeter 
equation J5j. This task will be accomplished in the future. 

Acknowledgements 

This article was supported by the DFG within the Sonderforschungsbereich 652 'Starke Korre- 
lationen und kollektive Phanomene im Strahlungsfeld: Coulombsysteme, Cluster und Partikel.' 
A.W. would like to thank the Center of Atomic and Molecular Technologies of Osaka University 
for its hospitality. C.F. acknowledges stimulating discussion with J. Vorberger and W.D. Kraeft 
and thanks R. Zimmermann for many helpful advice. 

References 

[1] H.A. Griem, Principles of Plasma Spectroscopy (Cambridge University Press, Cambridge, 
1997). 

[2] I.H. Hutchinson, Principles of Plasma Diagnostics (Cambridge University Press, Cambridge, 
1987). 

[3] see the proceedings of earlier workshops on the Physics of non-ideal Plasmas, e.g. Contrib. 
Plasma Phys. 39, 5-184 (1999); 41, 119-302 (2001); 43, 241-397 (2003). 

[4] M.D. Perry and G. Mourou, Science 264, 917 (1994); C.J. Joshi and P.B, Corkum, Physics 
Today 108, 36 (1995). 

[5] A.L. Fetter and J.D. Walecka, Quantum Theory of Many Particle Systems (McGraw-Hill, 
New York, 1971). 

[6] CD. Mahan, Many-Body Physics (Plenum, New York, 1990). 

[7] W.D. Kraeft, D. Kremp, W. Ebeling, and G. Ropke, Quantum Statistics of Charged Particle 
Systems (Plenum, New York, 1986). 

[8] J. Ward, Phys. Rev. 78, 182 (1950); Y. Takahashi, Nuovo Cimento 6, 370 (1957); T. Toyoda, 
Ann. Phys. (NY.) 173, 226 (1987). 

[9] G. Baym and L.P. Kadanoff, Phys. Rev. 124, 287 (1961). 

[10] In plasma physics, the use of V s instead of W is more common. 

[11] L. Hedin, Phys. Rev. 139, A796 (1965). 

[12] CD. Mahan, Comments Condens. Mater. Phys. 16, 333 (1994). 



11 



[13] F. Aryasetiawan and O. Gunnarsson, Rep. Prog. Phys. 61, 237 (1998); G. Onida, L. Reining, 
and A. Rubio, Rev. Mod. Phys. 74, 601 (2002). 

14] R.W. Godby, M. Schluter, and L.J. Sham, Phys. Rev. B 37, 10159 (1988). 

15] J.E. Northrup, M.S. Hybcrtson, and S.G. Louie, Phys. Rev. Lett. 59, 819 (1987). 

16] F. Aryasetiawan, Phys. Rev. B 46, 13051 (1992). 

17] E.L. Shirley, Phys. Rev. B 54, 7758 (1996). 

18] B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998). 

19] A. Schindlmayr and R.W. Godby, Phys. Rev. Lett. 80, 1702 (1998). 

20] Y. Takada, Phys. Rev. Lett. 87, 224602 (2001). 

21] C. Fortmann, H. Reinholz, A. Wierling, and G. Ropke, Proceedings of the 28th International 
Workshop on Consended Matter Theories, in print, arXiv: physics/0502051. 

22] L. Hcdin and S. Lundqvist, in Solid State Physics, edited by H. Ehrenreich, F. Scitz, and D. 
Turnbull, Vol. 23, p.l (Academic Press, New York, 1969). 

23] U. von Barth and B. Holm, Phys. Rev. B 54, 8411 (1996). 

24] A. Wierling and G. Ropke, Contrib. Plasma Phys. 38, 513 (1998). 

25] M. Hindgren and C.-O. Almbladh, Phys. Rev. B 56, 12832 (1997). 

26] H.J. de Groot, P.A. Bobbert, and W. van Hacringen, Phys. Rev. B 52, 11000 (1995). 

27] T.M. Rice, Ann. Phys. (NY.) 31, 100 (1965). 

28] D.J.W. Geldart and R. Taylor, Can. J. Phys. 48, 155; 48, 167 (1970). 

29] G.D. Mahan and B.E. Scrnclius, Phys. Rev. Lett. 62, 2718 (1989). 

30] F. Bechstedt, K. Tenelsen, B. Adolph, and R. Del Sole, Phys. Rev. Lett. 78, 1528 (1997). 

31] R.T.M. Ummcls, P.A. Bobbert, and W. von Hacringen, Phys. Rev. B 57, 11962 (1998). 

32] B. Tanatar and E. Dcmircl, Phys. Rev. B 62, 1787 (2000). 

33] LA. Nechaev and E.V. Chulkov, Phys. Rev. B 71, 115104 (2005). 

34] K. Shindo, J. Phys. Soc. Japan 29, 287 (1970). 

35] J. Vorberger, M. Schlanges, and W.D. Kracft, Phys. Rev. E 69, 046407 (2004). 

36] A. Wierling, Th. Millat, G. Ropke, R. Redmer, and H. Reinholz, Phys. Plasmas 8, 3810 
(2001). 

37] Faleev et al., Phys. Rev. B 74, 033101 (2006). 

38] T. Schmielau, G. Manzke, D. Tammc, and K. Hcnnebcrgcr, Phys. Stat. Sol. (b) 221, 215 
(2000). 

39] R. Schepe, T. Schmielau, D. Tamme, and K. Hcnnebcrgcr, Phys. Stat. Sol. (b) 206, 273 
(1998). 

40] G. Manzke, T. Schmielau, and K. Henneberger, Contrib. Plasma Phys. 41, 207 (2001). 

41] G. Onida, L. Reining, R.W. Godby, R. Del Sole, and Wanda Andreoni, Phys. Rev. Lett. 75, 
818 (1995). 



12 



[42] D. Blaschke, C. Roberts, and S. Schmidt, Phys. Lett. B 425, 232 (1998). 

[43] A. H611, CD. Roberts, and S.V. Wright, nucl-th/0604029. 

[44] R. Fchr and W.D. Kracft, Contrib. Plasma Phys. 35, 463 (1995); W.D. Kracft, R. Fchr, 
Contrib. Plasma Phys. 37, 173 (1997). 

[45] B.I. Lundqvist, Phys. Konden. Mater. 6, 193; 6, 206 (1967); 7, 117 (1968). 

[46] Ben Yu-Kuang Hu, Phys. Rev. B 47, 1687 (1993). 

[47] R. Fehr, PhD thesis (Greifswald, 1997). 

[48] P. Ziesche, accepted for publication in Ann. Phys. (Leipzig) (2006). 

[49] W.J. Karzas and R. Latter, Astrophys. J. Suppl. Scr. 6, 167 (1967). 

[50] C. Fortmann, R. Redmer, H. Rcinholz, G. Ropke, W. Rozmus, and A. Wierling, J. High 
Energy Dens. Phys. 2, 57 (2006). 

[51] H. Bethe and W. Heitlcr, Proc. Roy. Soc. London A 146, 83 (1934). 

[52] M.L. Ter-Mikaelyan, Dokl. Akad. Nauk SSSR 94, 1033 (1953). 



13 



<o[Ry] 




Figure 4: Contour plot of the self-consistent spectral function with simplified self-energy (Eq. Ij21jl 
at solar core conditions as function of frequency u> and momentum p. The spectral function is 
asymetrically broadened and shifted from the quasi-particle energy luqp = p 2 . At high momentum 
the spectral function converges into a sharp quasi-particle resonance at the free-particle energy. 



14 




Figure 5: Imaginary part of the self-energy as calculated from Eq. J2U- The first (dashed) and 
second (dotted) iteration as well as the convergent result (solid) are shown. Inset: Comparison 
of GW approximation and lowest order vertex correction (dashed-dotted), cf. Eq. (|22|) . The 
exchange diagram gives a correction of at most 20%. 
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co[Ry] 

Figure 6: Free-free absorption coefficient calculated from broadened electron propagators. The 
absorption coefficient is normalized to the Born approximation. At low frequencies, suppression 
of inverse bremsstrahlung is obtained. For high photon energies, the improved result converges 
to the Born approximation, while around lu = 10 Ry enhancement sets in. Also shown: a(u>) 
calculated with Gaussian spectral functions (dotted and dashed line) of different widths. Here, 
the convergence into the Born result is much slower, since the Gaussian spectral function does not 
yield the correct quasiparticle limit at high frequencies. Parameters: n — 7- 10 24 cm~ 3 , T = 100 Ry 
(solar core). 
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Figure 7: Absorption coefficient a{u) for solar core conditions as function of the photon energy 
u>. a(oj) is calculated from Eq. with broadened electron propagators using the self-energy 
obtained from Eq. 1|21|) (solid curve). Additionally, the first vertex correction is calculated, cf. 
Eq. I)22|l. The sum of both terms is given by the dashed curve. The vertex correction dominates 
the behaviour of the absorption coefficient. In particular, the enhancement at lu ~ 10 Ry present 
in the pure self-energy calculation vanishes completely. 
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